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Summary. — 

We improve the calculation of Nozawa, Itoh, and Kohyama (1998) to provide a 
formula for relativistic corrections to the thermal and kinematical Sunyaev-Zeldovich 
effects that is accurate to fourth order in 6 e = kBT e /m e c 2 , T e and m e being the 
electron temperature and electron mass, respectively. We also carry out a direct 
numerical integration of the Boltzmann collision term and confirm the excellent 
accuracy of the analytic formula. This formula will be useful for the analysis of 
the observational data of the forthcoming experiments of the kinematical Sunyaev- 
Zeldovich effect for clusters of galaxies. 

PACS 95.30.Cq - Elementary particle processes. 
PACS 95.30.Gv - Radiation mechanisms. 
PACS 98.65. Cw - Galaxy clusters. 
PACS 98.80-k - Cosmology. 



1. Introduction 

Compton scattering of the cosmic microwave background (CMB) radiation by hot 
intracluster gas — the Sunyaev-Zeldovich effect p]-[6] — provides a useful method for 
studies of cosmology (see recent excellent reviews: [7J, 0). Relativistic corrections to the 
Sunyaev-Zeldovich effect for cluster of galaxies have been extensively studied in recent 
years [9j-[23j. 

The present challenge to the observers of the Sunyaev-Zeldovich effect is the detection 
of the kinematical Sunyaev-Zeldovich effect, which is caused by the proper motion of 
galaxy clusters through CMB. By observing the kinematical Sunyaev-Zeldovich effect 
one can deduce the velocity of the galaxy cluster along the line of sight. Benson et al. 
[Mj have placed important upper limits to the velocities of the galaxy clusters along the 
line of sight with the method of the kinematical Sunyaev-Zeldovich effect. Within the 
next 5 years CMB experiments such as PLANCK, SPT, ACT, QUIET, APEX and AMI 
will perform deep searches for galaxy clusters with sensitivity limits at the level of 1-10 
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mJy. In the future CMB missions such as CMBPOL should reach sensitivities 20-100 
times better than those of PLANCK by using currently existing technology [35] . 

Under these current experimental circumstances it is worthwhile to obtain a preci- 
sion formula for the relativistic corrections to the kinematical Sunyaev-Zeldovich effect 
for clusters of galaxies. Nozawa, Itoh, & Kohyama [12] have calculated the relativistic 
corrections to the kinematical Sunyaev-Zeldovich effect for clusters of galaxies and have 
obtained a formula that is correct up to the order of {5Q\ where 9 e = kBT e /m e c 2 , T e and 
m e being the electron temperature and electron mass, respectively, and (3 = v/c is the 
magnitude of the peculiar velocity of the galaxy cluster with respect to CMB. Sazonov 
& Sunyaev [13] have calculated up to the (39 e term with a different method; their re- 
sult agreed with the (39 e term obtained by Nozawa, Itoh, & Kohyama [12] • Challinor & 
Lasenby [14j have carried out a further independent calculation and obtained the results 
that agreed perfectly with those of Nozawa, Itoh, & Kohyama [12] and Sazonov & Sun- 
yaev [13]. More recently Shimon & Rephaeli [21] carried out an independent calculation 
of the relativistic corrections to the thermal and kinematical Sunyaev-Zeldovich effects. 
However, their results for the kinematical Sunyaev-Zeldovich effect disagreed with those 
of Nozawa, Itoh, & Kohyama [12] . Sazonov & Sunyaev [15] ; and Challinor & Lasenby 
[14"] even at the level of the (39 e term. In this paper we will carry out the calculation 
up to the order of (39]; and obtain a precision formula for the relativistic corrections to 
the kinematical Sunyaev-Zeldovich effect for clusters of galaxies. We will also carry out 
a direct numerical integration of the Boltzmann collision term and confirm the excellent 
accuracy of our analytic formula. On the other hand, Shimon and Rephaeli's analytic 
formula shows considerable disagreement with the results of the direct numerical inte- 
gration of the Boltzmann collision term. Therefore we will conclude that Shimon and 
Rephaeli's analytic formula is in error. We will discuss on this in sec. 3. 

The present paper is organized as follows. In sec. 2 we briefly describe the method of 
the calculation and give the analytic solution. In sec. 3 we show the numerical results 
of the calculation. Concluding remarks will be given in sec. 4. In Appendix we describe 
our method of the numerical integration. 

2. Lorentz boosted Kompaneets equation 

In the present section we will follow Nozawa, Itoh, & Kohyama [12] and extend 
the Kompaneets equation to a system (the cluster of galaxies) moving with a peculiar 
velocity with respect to the CMB. We will derive higher-order terms than Nozawa, Itoh, 
& Kohyama [12]. We will formulate the kinetic equation for the photon distribution 
function using a relativistically covariant formalism. As a reference system, we choose 
the system that is fixed to the CMB. The z-axis is fixed to a line connecting the observer 
and the center of mass of the cluster of galaxies (CG). (We assume that the observer is 
fixed to the CMB frame.) We fix the positive direction of the z-axis as the direction of 
the propagation of a photon from the cluster to the observer. In this reference system, 
the center of mass of the cluster of galaxies is moving with a peculiar velocity (3{= v/c) 
with respect to the CMB. For simplicity, we choose the direction of the velocity in the 



In the CMB frame, the time evolution of the photon distribution function n(ui) is 
written as follows: 



x-z plane, i.e. (3 = ((3 x ,0,f3 z ). 
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(2) W 



(e 2 /Air) 2 X 5 i (p + k - p' - k 1 ) 



2ujlu'EE' 



(3) X^-(4 + ^ + w(I + l) 2 -W(i + l 

\K' K J \K K J \K K 

(4) k = -2(p • fc) = -2ujE (l - J-^cosa^ , 

(5) k' = 2(p • k') = UE (\ - ^cosa' 

In the above W is the transition probability corresponding to the Compton scattering. 
The four-momenta of the initial electron and photon arc p = (E,p) and k = (u>, 0, 0, k), 
respectively. The four-momenta of the final electron and photon are p' = (£" and 
k' = (uj',k'), respectively. The angles a and a' are the angles between p and k, and 
between p and k' , respectively. Throughout this paper, we use the natural unit Ti = c = 1 
unit, unless otherwise stated explicitly. Here we note that equation (1) assumes an 
isotropic photon distribution. The kinematic effects lead to anisotropy but to a good 
approximation equation (1) may be applied all the way through the cluster since the 
levels of anisotropy induced by the kinematic effect are small. 

The electron distribution functions in the initial and final states are Fermi-like in the 
CG frame. They are related to the electron distribution functions in the CMB frame as 
follows: 

(6) f(E) = .fc(E c ), 

(7) f(E') = f c (E' c ), 

(8) E C =J (e-P-p) , 

(9) E> c=1 (E> -p-i?) , 



(10) 



1 



where the suffix C denotes the CG frame. We caution the reader that the electron distri- 
bution function is anisotropic in the CMB frame for (3 ^ 0. By ignoring the degeneracy 
effects, we have the relativistic Maxwellian distribution for electrons with temperature 
T e as follows: 



fc(E c ) 
(11) 



e {(Ec-rn)-(p c -m)}/k B Te + j 
-{{E c -m)-(iJ.c-m)}/kBT e 



-1 



where (p,c — m) is the non-relativistic chemical potential of the electron measured in the 
CG frame. We now introduce the quantities 
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Substituting equations (6) - (13) into equation (1), we obtain 

= -2 



' , " U ' ) - * f ^ d 3 p'd 3 k'W fc(E c ) [{l+n(u/)}n(u) 



at 



(2tt) s 



(14) 



{ 1 + n(u) } n(w') e**rV-W) e ^-(k-~k') 



where k and k' are the unit vectors in the directions of k and k', respectively. Equation 
(14) is our basic equation. 

We now expand equation (14) in powers of Ax by assuming Ax <C 1. We obtain the 
Fokker-Planck expansion 



dn(u>) 
dt 



(15) 



where 
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(17) J ^ = ^j ^d 3 p'd 3 k'W fc(Ec) (Ax) k (l 



^ 7 /3-(fe-fc') 



Analytic integration of equations (16) and (17) can be done with the power series expan- 
sion approximation of the integrand in terms of the electron momentum p. In Itoh et 
al. [llj . the systematic analysis has been done in order to examine the accuracy of the 
power series expansion approximation. It has been found that the power series expansion 
approximation is sufficiently accurate for ksT e < 15keV by taking into account 0(9^,) 
corrections. 

In addition to 9 e , there is another parameter [3 in equations (16) and (17). For most 
of the cluster of galaxies, (3 <§; 1 is realized. For example, j3 « 1/300 for a typical 
value of the peculiar velocity w=l,000 kms -1 . Therefore it should be sufficient to expand 
equations (16) and (17) in powers of [3 and to retain up to 0((3 2 ) contributions. We 
assume the initial photon distribution of the CMB to be Planckian with a temperature 
T : 



(18) 
where 

(19) 



n (X) 
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Substituting the results of equations (16) and (17) into equation (15) and assuming 
T /T e -C 1, one obtains the following expression for the fractional distortion of the photon 
spectrum: 



An(X) 




rXe x 


n Q (X) 




e x - 1 
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rXe x 




e x - 1 




+ 


rXe x 




e x - 1 




+ 


rXe x 




e x - 1 



+ ^— y /3Pi(cos0 7 ) [ c + e e d + e 2 e c 2 + e 3 e c 3 + e 4 e d] 

tXp x 

(20) + -^—^(3 2 P 2 (cos9 1 ) [D + 6 e D 1 + 8 2 e D 2 + 8 3 e D 3 ] , 



(21) t = a T J d£N e , 

„ _ k B T e 



m e c 2 



(23) cos6» 7 



(22) e e = 

fa 

-1-0' 

(24) Pi(cos# 7 ) = cos# 7 , 

(25) P 2 (cos(9 7 ) = i (3cos 2 6» 7 - l) , 

where 7 is the angle between the directions of the peculiar velocity of the cluster (3 = v/c 
and the photon momentum k, which is chosen as the positive z-direction. The reader 
should remark that this sign convention for the positive z-direction is opposite to the 
ordinary one. Thus a cluster moving away from the observer has cos0 7 < 0. However, 
because of the positive sign in front of the Pi(cos# 7 ) term, one obtains An(X) < in 
this case, as one should. The coefficients are defined as follows: 

(26) Y = -4+X, 

(27) y 1 = -i 0+ £*-f* + ^ + *(-£ + |* 

Y 2 = - ±* + mx - ™X 2 + ™X* - ^X 4 + HjfB 
2 8 5 5 5 30 

^2 ( 434 658 - 242 ~ , 143 - ,\ ~ 4 / 44 187 - 

(28) +5 2 1 X X 2 -\ X 3 + S 4 \ 1 X 

1 ' V55 5 30 yV5 60 

15 2505- 7098-, 14253 -o 18594 ~ 4 

Y 3 = 1 X X 2 H X A X 4 

2 8 5 10 35 

12059 - 5 128 - 6 16 - 7 
+ ^4L> ~^l" X + 105 X 

§2 / 7098 14253 - 102267 - 2 156767 - 3 1216 ^ 64 ^, 
III 5 35 140 7 7 

18594 205003 - 1920 - 2 1024 - 3 

; 35 280 7 35 
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135 30375- 62391-, 614727 ~, 124389 

Y A = 1 X X 2 H X A X 4 

32 128 10 40 10 

355703 - 5 16568 - 6 7516 - 7 22 - 8 11 ~ 9 

80 21 105 7 210 

~ 2 / 62391 614727 - 1368279 - 2 4624139 - 3 157396 ~ 4 

30064 - 5 2717 - 6 2761 - 7 

7 7 210 

~ 4 / 124389 6046951 - 248520 - 2 481024 - , 15972 - , 

+ S 14 1 X X 2 H X 3 X 4 

V 10 160 7 35 7 

18689 - 5 

H X 5 

140 

^^f 70414 465992 - 11792 - 2 19778 - . 
+ 5 r^r + ^OS - ^"^^ + ^05~* 

(30, + + 

(31) B = ^o, 

(32) B 1 = ^Y + ^Yi , 

(33) B 2 = ^Y + ^Y 1 +Y 2 , 

(34) B 3 = -^Y + ^Y 1 + ^Y 2 + ^Y 3: 



(35) C = l. 



(36) CWO-^X + ^ + ^-S 2 , 
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'2 
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^ +s {-2o-lT x+ x) X ) + To s > 

75 21873- 49161~ 2 27519 -, 6684 ~ 4 
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.3917^5+64^ 
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+ 5 ^"SO 35~ X + ^5~ X ~ + ^5" X 
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^35 420 7 
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263259 - 5 4772 - 6 1336 - 7 11 - 8 
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30 \ 15 
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16364^ 5 _ 3764 ^ 6 101^ 
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4392 139094 - 3764 -, 6464 O 

h X X 2 H X 3 

5 105 7 105 j 



(43) + F(-*™ + ™*x) 1 
y ' V 315 315 J ' 

(44) l = Xcoth(^ , 



(45) 

sinh VT 



In the above, w is the photon frequency, To is the CMB temperature, <jt is the Thomson 
cross section, N e is the electron number density in the cluster frame, and the integral in 
equation (21) is over the photon path length in the cluster. 

Concerning the photon number conservation of equation (20), the following remarks 



8 



NOZAWA ET AL. 



should be noted: As for the thermal Sunyaev-Zeldovich effect terms, the photon number 
conservation is satisfied as discussed in Itoh, Kohyama & Nozawa [UJ. Namely, Yo, Y\, 
Y2, Y3, and I4 terms separately vanish by the integration J dXX 2 . Therefore the first 
and second lines vanish by the integration. On the other hand, the third and fourth 
lines are proportional to Pi(cos6> 7 ) and P2(cos0 7 ), respectively. Therefore the third and 
fourth lines vanish by the integration over the solid angle J c?f2 7 . 

Finally, we define the distortion of the spectral intensity as follows: 

(46) A/ = ^T~T~TFT 1 

e A - 1 n {X) 

Here we note that Itoh, Kohyama, & Nozawa [TT] have carried out the exact numerical 
integration of the Boltzmann equation for the thermal Sunyaev-Zeldovich effect and 
have compared the results with those obtained by the generalized Kompaneets equation, 
which correspond to the first line in equation (20). They have thereby confirmed that 
the analytic expansion result including up to the 6fYi term agrees with the exact result 
within the accuracy of about 1% for all the values of X for the electron temperature 
k B T e < 15keV (0 e < 0.03). 



3. Numerical results 



We now study the higher-order terms of the relativistic corrections to the kinematical 
Sunyaev-Zeldovich effect. In fig. 1 we have plotted the distortion of the spectral intensity 
AI/t as a function of X for fc s T e =10keV and /3 = /3 Z = 1/300. The thermal and 
kinematical Sunyaev-Zeldovich effects have been plotted. For the kinematical Sunyaev- 
Zeldovich effect, the absolute value of the effect is shown. For a receding cluster one 
has AI < 0, whereas for an approaching cluster one has AJ > 0. We have shown the 
decomposition of the absolute value of the kinematical Sunyaev-Zeldovich contribution 
in fig. 2. The dotted curve is the leading order contribution originally derived by Sunyaev 
& Zeldovich [4] , [5] . The dashed curve is the result of Nozawa, Itoh, & Kohyama [12] , 
which includes up to O(f30 2 ) terms. The dot-dashed curve is the present result that 
includes up to 0((38l) terms. The solid curve is the present result that includes up to 
0(/36f) terms. The triple-dot-dashed curve is the contribution of the 0({3 2 ) terms, where 
a factor 10 was multiplied in order to make the curve visible in this figure. It is clear 
from fig. 2 that the convergence of the power series expansion is extremely good at this 
temperature for all X values. It is also clear from fig. 2 that the 0(/3 2 ) correction is very 
small and it can be safely neglected. 

It is well known that the kinematical Sunyaev-Zeldovich effect becomes extremely 
important in the crossover frequency region, where the thermal Sunyaev-Zeldovich effect 
vanishes. Therefore the accurate position of the crossover frequency Xq is extremely 
important. In Itoh, Kohyama, & Nozawa [UJ, the following fitting formula for Xq has 
been obtained. 



(47) Xq = 3.830 (1 + 1.16740 e - 0.8533^) . 

The errors of this fitting function were less than 1 x 10~ 3 for < fcsT e < 50keV. 

Substituting equation (47) into the third line of equation (20), we have calculated the 
kinematical Sunyaev-Zeldovich effect at X = Xq as a function of ksT e for (3 — /3 z =l/300. 
The result has been plotted in fig. 3. The dotted curve is the result of Nozawa, Itoh, 
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Fig. 1. - Spectral intensity distortion AJ/r as a function of X for ksT^ — lOkeV, /3 = /3 g = 1/300. 
The thermal and kinematical Sunyavev-Zeldovich effects have been plotted. For the kinematical 
Sunyaev-Zeldovich effect, the absolute value of the effect is shown. 
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Fig. 2. - The decomposition of the absolute value of the kinematical Sunyaev-Zeldovich contri- 
bution for fcsT e = lOkeV, /3 = /3 z =l/300. The dotted curve is the leading order contribution 
originally derived by Sunyaev & Zeldovich [3], |5]. The dashed curve is the result of Nozawa, 
Itoh, & Kohyama |12j . which includes up to O(PO^) terms. The dot-dashed curve is the present 
result that includes up to O(f30'i) terms. The solid curve is the present result that includes up 
to 0(f3@i) terms. The triple-dot-dashed curve is the contribution of the 0((3 2 ) terms, where a 
factor 10 was multiplied in order to make the curve visible in this figure. 
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Fig. 3. - The absolute value of the kinematical Sunyaev-Zeldovich effect at X = Xo as a function 
of fcsTe for (3 — f3 z = 1/300. The dotted curve is the result of Nozawa, Itoh, & Kohyama [12] , 
which includes up to O(f30'l) terms. The dashed curve is the present result that includes up to 
0(f38i) terms. The dot-dashed curve is the present result that includes up to 0([38f) terms. 
The solid curve is the present result of the direct numerical integration. 



& Kohyama [12], which includes up to 0((38 2 ) terms. The dashed curve is the present 
result that includes up to O(f30l) terms. The dot-dashed curve is the present result that 
includes up to 0(08*) terms. The solid curve is the present result of the direct numerical 
integration. It is clear from fig. 3 that the convergence of the power series expansion is 
extremely good for the temperature ksT e < 20keV. 

We have derived an analytic expression for the Sunyaev-Zeldovich effect in the Rayleigh- 
Jeans limit where X — > 0. The result is 



An(X) 

MX) 



(48) 



-2r( 



-2t(3 2 



l- l lt 
10 
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— i 
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+ t/3J>i(cos(9 7 ) 
+ t/? 2 P 2 (cos6> 7 



40 e 
23. 
' + 20* 

5 

1 4 
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15 5 



1989 
280 ' 
2539. 
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13 -a 

T 0e 



7 



14403 ( 
640 



1689^3 
140 e 

2 341 

e ~ IT* 



7281 

( 

140 



We have also carried out a direct numerical integration of the Boltzmann collision 
term on the right-hand-side of equation (14). The method of the numerical integration 
is presented in Appendix. Here we will show an example of our numerical calculations. 
We show the results of the following case. We have considered the case that # e =0.02 
and the cluster approaches the observer (cos# 7 =l) with peculiar velocity u=1000 kms -1 . 
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, = 0.02, (3 = (3 =1 /300 



-2 - 



O 



-4 



-6 



15 



Fig. 4. - The comparison of 5 P rcscnt (solid curve) and 5sr (dashed curve) as defined by equations 
(49) and (50), respectively. 



The numerical result for the kinematical Sunyaev-Zeldovich effect has been obtained by 
subtracting the numerically-obtained purely thermal effect (/?=0) from the numerically- 
obtained total (thermal plus kinematical) effect. Then we have compared this with 
our formula for the kinematical Sunyaev-Zeldovich effect (the sum of the second, third, 
and fourth lines of equation (20)) and also with Shimon and Rephaeli's formula for the 
kinematical Sunyaev-Zeldovich effect. 
We define the following quantity 



(49) 
(50) 



^present — 
SSR = 



(^-Oprcscnt analytic (^-Opresent numerical 



(^-Oprcscnt numerical 
(A/)sR analytic ~ (^-Oprcscnt numerical 



(^-Oprcscnt numerical 



The results are shown in fig. 4. It is readily seen that our analytic formula for the kine- 
matical Sunyaev-Zeldovich effect agrees with the numerical results for small values of 
X with accuracy better than 10~ 4 . On the other hand, Shimon and Rephaeli's formula 
shows about 2% discrepancy with the numerical results for small values of X. Therefore 
we have confirmed the excellent accuracy of our analytic formula for the kinematical 
Sunyaev-Zeldovich effect. We conclude that Shimon and Rephaeli's formula for the kine- 
matical Sunyaev-Zeldovich effect is in error. 

A discussion on the discrepancy of Shimon and Rephaeli's results from ours is in 
order. All of our calculations have been carried out consistently in the CMB frame. 
The electron distribution function f(E) in the CMB frame in equation (1) is equal to 
fc(Ec), the electron distribution function in the CG frame (equation (6)). The electron 
energy in the CG frame Eq is expressed in terms of the electron energy and momentum 
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in the CMB frame (equation (8)). The integration is carried out consistently over the 
momentum space in the CMB frame (equation (1)). However, Shimon and Rephaeli 
appear to have erroneously used the momentum space in the CG frame instead. This is 
the likely cause of the discrepancy of Shimon and Rephaeli's results from ours. Shimon 
and Rephaeli correctly reproduce the Co term in equation (20), and their discrepancy 
starts at the 9 e C\ term. Since the contribution of the 9 e C\ term with respect to the 
contribution of the Co term for the case of 6*e=0.02 is about 2% for small values of X, 
this explains the results in fig. 4. 

4. Concluding remarks 

We have extended the previous calculation of Nozawa, Itoh, & Kohyama (1998) on 
the relativistic corrections to the kincmatical Sunyaev-Zeldovich effect for clusters of 
galaxies to the order of (39* where 9 e = kBT e /m e c 2 . We have also carried out a di- 
rect numerical integration of the Boltzmann collision term and numerically calculated 
the kinematical Sunyaev-Zeldovich effect, thereby confirming the excellent accuracy of 
our analytic formula for the kinematical Sunyaev-Zeldovich effect. This precision for- 
mula should be very useful for the analysis of the observational data of the kinematical 
Sunyaev-Zeldovich effect for clusters of galaxies, which will be very likely to be obtained 
through the forthcoming CMB experiments. We have found that Shimon and Rephaeli's 
formula for the kinematical Suyaev-Zeldovich effect contains about 2% error for the case 
6* e =0.02, (3 = f3 z = 1/300 by comparing with the results obtained by the numerical inte- 
gration of the Boltzmann collision term. 

This is the first time that the accuracy of the analytic formula for the kincmatical 
Sunyaev-Zeldovich effect has been quantitatively assessed by comparing with the results 
obtained by the numerical integration of the Boltzmann collision term. It is gratifying 
to find that the excellent accuracy of the analytic formula for the kinematical Sunyaev- 
Zeldovich effect has been proved by comparing with the numerical results. 
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Appendix 



In this appendix we describe our method of the numerical integration of the Boltz- 
mann collision term on the right-hand-side of equation (14). Let us first start with the 
electron distribution function in the CG frame, i.e. equation (11). 



fc(Ec) = 



(A.l) 



e {(Ec-m)-( MC -m)}/k B T e + ^ 
-{(Ec-m)-(pc-m)}/k B T s 



The electron number density N e is related to the electron distribution function as follows. 



Nf, = 2 



d 3 p c 



fc (E c ) 



(A.2) 



o W 
1 r°° 

= J_ e (»c-m)/k B T e / dpc ^ e - 

n Jo 



K c /k B T c 



where Kq = Ec — m. This equation is further simplifcd. 



N e = e (^-m)/feBT e ™ ^ r du (1 + (2 + ^ e -u 

n Jo 

Therefore the chemical potential term in fc(Ec) is expressed by N e as follows. 

(A.3) e {nc-m)/k B T e = Ne f_Q~i g (Q e \ ; 

" poo "I — 1 

(A.4) g(0 e ) = / dx(l + 9 e x) ^Jx (2 + 9 e x) e~ x 

Jo 

The equation (14) is rewritten with equation (2) as follows. 

dn(uj) 



dt 



= -2 



d 3 p d 3 p' d 3 k'X sit , 

——(5 (p + k-p ~k)f c (E c ) 



4tt J J (2tt) 3 2lulu'EE' 
x |n(w) [1 + n (lu')] - n (lu') [1 + n(u)] e ^(^-0-h') e ^0-(k-k') j 



= -2 



d a p 1 
(2tt) 3 £^ 



J dVL k , J dw'J 5 (p + k- k'f - to 2 ] f c {E c ) X 



(A.5) 



{n(w)[l + n (lu')} - n (lu') [1 + n(Lu)]e Ax ^^' k ') e x ^-' k ') } 



Equation (A.5) is further simplified as follows. 



dn(u) 



h) Iwphl ^ 7u~ E + cu ( 1 - cos 9 Y ) p cos % fG {Ec) * 



(A.6) x {n(co) [1 + n (lu')] - n (lu') [1 + n(w)] e A *7(i-/3-fc') e *7£-(fc-fc') } 
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and 

(A 7) J = u(E~pcos9 p ) 

{ '' E + u(l - cosfly) -pcos6' p ' 

where 9 P , 9' p and 9 1 ' are the angles between p and k, p and A;', fc and fc', respectively. In 
equation (A. 6) the chemical potential term of fc(Ec) is replaced by equation (A. 3). 

Finally we perform the five-dimensional integration of equation (A. 6). We use the 
Gaussian quadrature integration method, which is commonly used for the numerical 
integrations. The method is very efficient as well as accurate. We have checked the 
convergence of the integral carefully by varying the number of quadrature points. The 
accuracy of the integral is 10~ 6 for the case of 70 quadrature points for each integral 
variable. The numerical integration itself is straightforward, although it is rather time 
consuming. It takes about 20 hours (CPU time) for a standard Alpha Server Workstation 
in order to complete the calculation of the spectral intensity AI for fixed ((3, 9 e ) values 
and for < X < 20 with 100 points. 
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